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We present the results for the hadron spectrum calculated on 400 configurations using point source, wall source 
and 8-cubic sources, in quenched QCD with Wilson fermions at = 6.0 and K = 0.155 on a 24 3 x 54 lattice. The 
results for the ground state masses obtained with three types of quark sources agree well with each other. Masses 
of the first excited states appear consistent with experimental values within large errors. 
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Recent high statistics simulations of quenched 
QCD have produced high quality data for the 
masses of light hadrons[0-||. However, mass re- 
sults obtained with various types of quark sources 
often do not agree with each other HH^. Masses 
of the first excited states also often appear unex- 
pectedly heavy dJ|J^l. 

In order to resolve these problems, the QCD- 
PAX Collaboration has carried out a high statis- 
tics quenched calculation with Wilson quarks at 
0= 6.0 on a 24 3 x 54 lattice. Taking the last con- 
figuration of our previous work|l| as a starting 
configuration, we have generated new 400 config- 
urations which are separated by 1000 heat bath 
sweeps. We choose K = 0.155: This K corre- 
sponds approximately to the strange quark mass, 
because the rn T /m p appears to be 0.7. 

We use three types of quark sources: point 
source, wall source and 8-origin cubic sources of 
a size 7 3 (with a random Z(3) element on each 
cube). We call the third source multi-cube source. 
Hereafter we abbreviate them to p-, w- and m- 
sources, respectively. Point sinks are used for 
all sources. Configurations are fixed to Coulomb 
gauge before quark propagators are solved for the 
w- and the m- sources. Although we have 600 
samples for the p-source in total, we report the 
results obtained by using the new 400 configura- 
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tions, in order to compare results with three types 
of sources on an equal footing. 

Masses of the ground states are determined 
by global fits to data for [t m i n ,t max ] for which 
the x 2 /dof of the fits are less than 1.5. We 
take account of time correlations in the fits. Er- 
rors are estimated by single elimination jack-knife 
procedure. We have checked that errors do not 
change significantly even if we change the bin 
size in the jack-knife error analyses. In addition, 
we have estimated errors which depend on the 
choice of fitting ranges. Here we define the fitting- 
dependent upper (lower) bound by the maximum 
(minimum) of fit results to data for any region 
[to,h] C [t m in,t ma x\ (h - *o > 2). Detail of our 
fitting procedure will be reported elsewhere. 

Results for the ground state masses are repro- 
duced in Table [l] and are summarized as follows. 
1) Long plateau in effective masses can be seen 
for all particles and all sources. Length of the 
plateau (t rnax - t mm + 1) is typically 15 for tt, 
8 for p, 6 for nucleon and 8 for A. 2) Statistical 
and fitting-dependent errors are typically of order 
2 %. 3) Mass results obtained with three types 
of quark sources coincide within 0.3 % for tt, 1 
% for p, 4 % for nucleon and 2 % for A. They 
agree within errors. (See also Figs, [j] and |[) It 
may be worth while emphasizing that the A mass 
obtained with the p-source agrees with that with 
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Table 1 

Results for the ground state masses. In parentheses are errors estimated by the jack-knife method. Errors 
given in the form ^ arc for the fitting dependent upper/lower bound. 
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Multi-Cube Source 
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0.2960( 8) 
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0.4201(29) 
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24 


0.6403(50) 
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23 


0.7054(95) 
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trnax 


mass(er) 
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- 27 


0.2964( 6) 


13 


- 20 


0.4228(19) 


13 


- 18 


0.6307(39) 


11 


- 18 


0.7008(57) 
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mass(er) 
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+ 9 

- 20 
+ 25 

- 36 
+ 
-161 
+ 9 
-171 



12 
5 
8 
15 
44 
10 
14 
59 



12 - 27 
11 - 20 
14-18 
14-24 



0.2969( 6) t J 
0.4249(19) ± 33 
0.6208(72) ± l° 3 

+ 13 

-345 



0.7128(191) 



the w-source with small errors. Note that the er- 
ror for the A mass is the largest for the m-source. 
The agreement of the results obtained with vari- 
ous types of sources implies that they have little 
contamination from excited states. (Sec also the 
contributions by Ape group|^| and R. GuptaQ.) 
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Figure 1. p meson masses obtained with vari- 
ous types of quark sources. Two error bars are 
drawn for our new results (filled circles). On 
the left are for errors estimated by the jack-knife 
method. The right ones are for fitting depen- 
dent upper/lower bounds. Other data are from 
LANL group§ (diamonds), Ape Collab.|] (tri- 
angle) and our previous workjl| (open circles). 



The results for hadron masses are consistent 
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Figure 2. The same as fig. [l] but for A. 



with our previous ones. In ref.§, we have made 
two different fits (plateau fit and pre-plateau fit) 
to our old p meson propagator and obtained m p = 
0.4169(48) and 0.4280(33) for each fit, respec- 
tively. We have taken the plateau fit to be cor- 
rect. Our new results are consistent with the 
value obtained by the plateau fit. (See fig. [[].) For 
baryons, new results are consistent with the old 
ones which are obtained by two-mass fits; m^— 
0.635(11) and m A = 0.711(12). Note that in our 
previous work the errors of baryon propagators 
for the p-source are large at high time slices and 
consequently the one mass fit gave a larger value 
than the two mass fit. We have taken the re- 
sult by the two mass fit as correct. In this work 
the errors have been reduced significantly due to 



3 



high statistics. Therefore the results of one mass 
fits are smaller than our previous ones (mjv= 
0.645(6) and toa= 0.728(7)) and are consistent 
with our previous results of two mass fits. 

Results for the mass ratios m^/mp and 
mA/ffl p are consistent with values estimated by 
phenomenological mass formulae JTo| : mjv/m /9 = 
1.524 ±47, 1.492 ±^, 1.461 ±}g and m A /m p = 

1.679 ±51,1.658 ±",1.678 ±gg for the p-, w- 
and m- sources, respectively. We have estimated 
two kinds of errors. One is the error obtained 
by the jack-knife method for the mass ratios. 
The other is the fitting dependent error estimated 
by dividing the fitting dependent upper (lower) 
bound for the numerator by the fitting depen- 
dent lower (upper) bound for the denominator. 
We have quoted the larger one. 

Our next question is whether the masses of the 
first excited states appear consistent with exper- 
iment. We find that for each particle there exist 
two mass fits with reasonable \ 2 ' s which give the 
ground state masses consistent with one mass fits. 
For these fits, almost all the results for the first 
excited state masses are consistent with experi- 
ment within large one standard deviation. 

Result for the p meson is summarized as fol- 
lows. For the p-source, the fit with t m in = 7 
gives to p = 0.4184(34) and m' p = 0.836(63) with 
X 2 /dof ~ 1.7. For the w-source, all the fits with 
train — 4^7 give x 2 /dof — 1 and m' p ~ 0.7. We 
quote the result from the fit with t min — 6: m p = 
0.4240(30), toJ=0.72(11). The results for m' p are 
consistent with the experimental value m' p ~ 0.7. 
(This value is obtained assuming that the K used 
corresponds exactly to the mass of the strange 
quark and using the experimental values for <\> me- 
son masses; 0(1020), 0(1680).) However, for the 
m-source, there are only fits which are marginally 
consistent with m' p ~ 0.7 with large errors, for ex- 
ample, 0.96(33) with t min — 8. 

For nucleon, the fit with t m . m = 8 (7, 6) for the 
p- (w-, m-) source gives, with % 2 /dof =1.1 (0.8, 
1.2), m N = 0.6302(77) (0.6346(79), 0.6253(128)) 
which is consistent with the one-mass fit result. 
The results for the m' N agree with each other 
with large errors: m' N = 0.99(13), 0.93(13) and 
1.01(7), respectively. We estimate m' N — 0.8 ~ 



1.1 from these fits. This value is smaller than that 
reported by Ape group (1.38(2))||. We estimate 
the experimental value of the m' N to be 0.84 as- 
suming that m' N — tojv at the strange quark mass 
is the same as that at u/d quark mass (500 MeV). 

Two mass fits are more delicate for A. With 
the constraint of reasonable x 2 , only the fit with 
tmin — 9 for the p-source gives the result for the 
m/\ which is consistent with that obtained by one- 
mass fit. This fit gives m\ — 0.695(14) and m' A = 
1.02(11) with x 2 /dof=l.l. Experimental value 
for the m' A is 0.9 (from fi(2250)"). 

The numerical calculations have been per- 
formed with QCDPAX, a parallel computer 
developed at University of Tsukuba. This 
project is supported by the Grant-in-Aid of 
Ministry of Education, Science and Culture 
(No.62060001,No.02402003 and No.04NP0601) 
and University of Tsukuba Project Research in 
1993. One of the authors (T.Y.) would like to 
thank Yamada Science Foundation for financial 
support. 
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